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ABSTRACT 



We present a detailed theoretical analysis, using correlation functions, of the coherence properties of the 
output from a frequency shifted feedback (FSF) laser seeded simultaneously by an external seed laser and 
by spontaneous emission (SE). We show that the output of a FSF laser is a cyclostationary process, for 
which the second-order correlation function is not stationary, but periodic. However, a period-averaged 
correlation function can be used to analyze the optical spectrum. From the fourth-order correlation func- 
tion of the output of a Michelson interferometer we obtain the essential characteristics of the radio-fre- 
quency (RF) spectrum, needed for describing the use of the FSF laser for optical-ranging metrology. We 
show that, even for a FSF laser seeded by SE, the RF spectrum comprises a sequence of doublets, whose 
separation gives directly a measure of the length difference between the interferometer arms. This dou- 
blet structure is a result of the correlation of interference terms of individual components of the cyclo- 
stationary stochastic process. It is not seen in the optical spectrum of the FSF laser but is observable in 
the RF spectrum. We analyze the competition between SE and continuous wave (CW) seeding to obtain 
an analytical expression for the ratio of power in the discrete CW signal to the background continuum 
spectrum from SE. We show that, unlike mode competition in conventional lasers, where there occurs 
exponential selectivity, here there is a balance between the two fields; the power in the fields is related 
linearly, rather than exponentially, to the control parameters. 

© 2008 Elsevier B.V. All rights reserved. 



1. Introduction 

The remarkable properties of frequency shifted feedback (FSF) 
lasers have drawn attention as broad-band light sources [1-4] 
and, more recently, as tools for metrology [5-9]. In essence, an 
FSF laser comprises a closed optical path, of total optical length L 
(a linear cavity or a ring), in which light undergoes not only the 
usual gain of energy from excited states but also, during each cav- 
ity round trip, a discrete-frequency shift A. This differs from the ax- 
ial mode spacing of the cavity, the free-spectral range (FSR). which 
is 2ttc/L. The operation of FSF lasers has been considered by several 
research groups [ 1 0-20] and it has been used for a number of prac- 
tical applications [21-26]. 

Two complementary viewpoints offer means for understanding 
the output characteristics, and uses, for an FSF laser. In one. the 
output field is regarded as a comb of chirped frequencies, sepa- 
rated in frequency at any instant by the free-spectral range 
2nc/L (the moving-comb model) [5.6,14-1 6.27]. The other approach 
considers the output growing from a discrete-frequency seed as a 
stationary set of frequencies starting from the seed frequency 
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and separated from each other by A [19.20,23.28.29]. Spontaneous 
emission (SE) provides a continuous distribution of such seed 
frequencies. 

For applications to metrology, specifically the determination of 
distances by optical ranging, it is important to understand the 
coherence characteristics of the output from an FSF laser. In the ab- 
sence of any external seed, this output originates from a distribu- 
tion of spontaneous-emission (SE) events within the gain 
medium. These are uncorrelated and intrinsically incoherent. We 
here examine, by following a sequence of cyclic passages within 
the FSF cavity, how this radiation acquires coherence properties, 
as evidenced in correlation functions. 

We discuss here the operation of a FSF laser whose field origi- 
nates either entirely from spontaneous emission or else from the 
joint action of spontaneous emission and a continuous-wave 
(CW) monochromatic laser seed. We pay particular attention to 
the RF spectrum of the output of a Michelson interferometer driven 
by a FSF laser. As we will note, this signal underlies metrology 
applications. 

Earlier we have presented the basic theory of the operation of a 
FSF laser [19]. That work treated the SE as a seed field of very large 
bandwidth. It also described the effect of an external CW laser seed 
field, one whose phase could be modulated. However, it did not 
discuss the competition between these two sources. Here we rem- 
edy these shortcomings by treating the two sources as coexisting. 
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All applications of FSF lasers to distance measurements {optical 
ranging) rely on measurements of the RF spectrum of the output of 
a Michelson interferometer with FSF-Iaser input. Hitherto there 
has been no quantitative ab initio theory of the properties of this 
spectrum. Such a description must originate with intensity correla- 
tion functions, as presented here. We will show how the fre- 
quency-shifting element produces a field which, when passed 
through a Michelson interferometer, produces spectral signals 
from which the interferometer-arm distance can be deduced. Spe- 
cifically, the information is contained in sets of frequency doublets. 
We present here a first-principles derivation of the "beat frequen- 
cies" that have been discussed earlier based on the moving-comb 
model of FSF operation. We show that such signals occur even 
when the FSF laser is seeded entirely by spontaneous emission, 
and that there is no need to invoke a moving comb model to ex- 
plain this structure. The interferences from such SE sources are 
not seen in the optical spectrum of the FSF laser; they are observa- 
ble only in the RF signals from the interferometer. 

Our earlier work [19] considered the FSF output field from a 
specific seed field. We used these results to model the output from 
an FSF laser seeded either solely by SE or exclusively by an exter- 
nally controlled laser field. As shown in (7), the latter option offers 
significant improvement in signal-to-noise ratio over SE seeding. 
The present paper treats FSF-Iaser operation when an external seed 
field is accompanied by the inevitable spontaneous emission that 
provides an unavoidable stochastic background. In normal laser 
operation the competition between modes leads to exponential 
growth of favored modes at the expense of those with less gain. 
We shall show that in the FSF laser the competition between SE 
and external seed does not act in this way. Instead, there occurs 
a balance between the two fields that is related linearly, rather 
than exponentially, to the control parameters. Although it would 
be desirable to suppress the amplified spontaneous emission 
(ASE) component of the FSF-Iaser output, thereby improving the 
signal-to-noise ratio, we shall show that this is not possible. This 
behavior differs, perhaps unexpectedly, from that of a normal laser, 
where suitable choice of gain characteristics can place all of the 
output into a favored mode. 

2. Basic model 

This article draws upon two recent papers [7,19] that present, in 
some detail, an idealized model for an FSF laser. We here first re- 
view the basic principles of that work and then extend it by exam- 
ining correlation functions, thereby allowing analysis of the 
coherence and spectral properties of the FSF laser. 

2.1. The FSF-Iaser layout 

Fig. 1 presents a schematic drawing of the essentials of a FSF la- 
ser in a ring configuration. Light propagates along a closed optical 
path that passes through a gain medium (G). a frequency shifter 
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[usually accomplished by an acousto-optic modulator, (AOM)] 
and a filter (F). In G a combination of stimulated emission and 
spontaneous emission add energy. The AOM adds a fixed incre- 
ment A to each frequency component. A number of elements with- 
in the overall optical path have frequency dependent losses that 
we here represent as a single spectral filter. F . The output emerges 
through the zero diffraction order of the AOM. The AOM also serves 
to couple any external seed laser beam into the cavity. 

2.2. The field description 

The filtering elements (F). which remove energy, and the gain 
medium (G). which replenishes this energy, together provide an 
effective gain (the growth minus the loss) distributed over a range 
of frequencies. We take the center of this range to be aj ft regarding 
this as a carrier frequency. We consider the propagation of the 
electric field along a one-dimensional path, coordinate z, within 
the ring. The round-trip time in the cavity is T r = L/c. 

We denote by E(r) the complex-valued electric field at position 
z = 0, taken to be the entrance surface of the AOM. We express this 
field in terms of a com pi ex- valued envelope (t) and a carrier: 

E(t) = 4(t)exp(-ia> f t). (1) 

The intensity of the FSF laser at z = 0, is 

/(0=^(t)| 2 . (2) 

The AOM diffracts into first order a fraction of the beam inten- 
sity, denoted Typically s* 2 % 0.9. The frequency of this fraction 
is shifted, in passing through the AOM, by A. This portion continues 
to circulate within the cavity. The unshifted fraction (1 - M 2 ), usu- 
ally small, remains in zeroth order and forms the output field. 

We consider an arbitrary but fixed reference time T and a time 
window that endures for one round-trip time prior to 7. Within 
this interval, T - T r < t < 7\ we introduce a sliding windowed Fou- 
rier transform (SWFT) of the envelope, writing 



(3) 



(4) 



& (f) = / dxD&{m, T) exp(-itnt), for T - x r < t < T, 

J -ao 

where 

1 f T 

4(m, T)= — J^ dtS(t) exp(+itnt). 

Here and in the following we denote by w = co - io f a frequency 
offset from the carrier co f . The SWFT field 6 (w, T) for a fixed refer- 
ence time T provides the essential tool with which we express the 
operation of the FSF and the output field [ 1 9], From it we construct 
the field envelope &(t) for all time, and deduce the output charac- 
teristics of the FSF for any given experimental conditions. 

2.3. The field amplitude equations 

As in our earlier work [19], we derive the basic equation for the 
laser field by considering its history during the round trip prior to 
time T. It began this circuit, starting at time t = T - r r , with fre- 
quency a) = co/ + w - A. It then underwent a frequency shift A as 
it passed through the AOM. It changed amplitude in passing 
through the filter and gain medium. It also acquired a phase incre- 
ment (a} f + w)t t as it propagated around the loop. To describe 
these changes it proves useful to introduce a two-dimensional 
time-frequency spaced = (tu, T) in which to present the equations 
[19]. Successive circulations around the ring involve a change of 
the two-dimensional vector X by the increment X 0 = (A. r r ). The 
field change after one round trip is expressible as 



Fig. 1. Symbolic diagram of a ring cavity showing gain G. spectral filter F <md 
frequency-shifter AOM as well as seed and output. 



S (X) = if (X - Xo) exp [C(X)J + {(X) + e(X). 



(5) 
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Here the effect of gain and loss appears through the propagator 
argument C(#). described in detail below, and the terms {(X) and 
£(X) represent, in Fourier space, additions to the field during the 
time interval r r . The first of these. £(*), expresses the field incre- 
ment supplied by spontaneous emission. It is the SWFT of 4s P (t). 
the amplitude of the stochastic field created by spontaneous-emis- 
sion events. 



(6) 



i(X) = i(^)=^J T \ *&p(0exp(+tot). 

The second contribution to the field. £{#). describes the external 
seed field, if present. We take this to have a carrier frequency oj s 
offset from the central filter frequency a) f by tn„ 

to, = to/ + nj s . (7) 

We are interested in seeding by monochromatic light; we write its 
SWFT as 



e(X) = I dt exp [-im s t + imt], 



(8) 



where e s is the amplitude of the seeding laser field within the 
cavity. 

We write the propagator argument C(#) as the difference be- 
tween gain g{co) and losses f{co), modified by the phase 
((Of + w)x, acquired during the round trip. We approximate the ef- 
fect of losses by means of a filter function, taken to be a quadratic 
centered at co = co f with characteristic width r f , 

rm-h +K^)'- « 

The frequency-independent term f m describes fixed losses, such as 
those from mirror reflectivities. 

We will regard the spectral filter as having narrower bandwidth 
than the broader bandwidth of the gain medium; it fixes the range 
of frequencies we need to consider. In this approximation we ne- 
glect the frequency dependence of the gain and approximate it as 

g (co)= gsat = -$—. (10) 

I + «/ 'sat y 

Here g sat is the saturated frequency-independent gain. g 0 is the 
unsaturated gain. 7 is the laser intensity averaged over one round 
trip and l SQt is the saturation intensity. The averaged intensity, and 
hence the saturated gain, must be evaluated from appropriate equa- 
tions describing the FSF laser, see Section 5.1 below. Combining 
these three contributions, we write the propagator argument as 



G(X) = 0(07, T) = i((o f + w)x r + g Ml - f m - ^ (w/r f ) 2 



This depends on frequency o through the detuning uj~co-co f . 

The basic parameters under the control of the experimenter are 
those of the cavity r r . J. those of the filter a) f . r f ,f m . those of the 
gain g 0 . I sat , and those of the seed. £„ a) s . The mean intensity 7. 
and the concomitant saturated gain g sat must be determined such 
that the equations for the field, with predicted output intensity, 
are self consistent. 

2.4. The output intensity 

An important characteristic of the FSF-laser output is the inten- 
sity averaged over a round trip, 



7(7) si/ dt/(f). 

T r Jt-t, 



(12) 



We define the spectral intensity l(xnj) by writing this average as a 
frequency distribution, 



/(T)= / dml(mj). (13) 

Using Eq. (3) we write the intensity as 

/(f) = ^ J j dxn x dxo 2 exp[-i(G7, - nj 2 )r)] 

xS(m } J)S(w 2 J)\ (14) 

From this integral we extract an expression for the spectral inten- 
sity of Eq. (13) as 



1 - exp -107'Tr 
x ^ ^ exp iaf 71 



(15) 



The fields described here, and the corresponding output inten- 
sity, involve uncontrollable events, the stochastic processes that 
describe spontaneous emission. Observations inevitably introduce 
averages over stochastic realizations of the average intensity dur- 
ing a round trip. Thus we are interested in stochastic averages 
< ) that give the spectrum 

J(m) = (l(m,T)), (16) 
and the mean intensity 

J=(/(T)). (17) 

2.5. The general solution for the output field 

The FSF laser based on the growth of spontaneous emission acts 
as a regenerative amplifier of spontaneous emission, as described 
by basic Eq. (5). The solution to this equation, as discussed earlier 
[19], is expressible as the contribution of two parts, 



a field that grows from the spontaneous emission. 
<M*) - «*) + £«* - n*o)«p[JU(*)l. 

n=l 

and a field that grows from the external seed, 
<WX) = e(X) + Je(X - nX 0 ) exp[A„(X)|. 
Here 



(11) a„(X) = jTC(X-lXo). 



(18) 



(19) 



(20) 



(21) 



(=0 



These formulas express the field as the most recent added incre- 
ments from seed and spontaneous emission, together with the suc- 
cession of fields from earlier passes, each altered by the 
appropriate propagator exp(;. n (X)j. Our earlier work considered 
each of these separately. Here we consider the operation of the la- 
ser when both are present. 

3. Seeding by spontaneous emission. The optical spectrum 

Earlier work on the theory of the FSF laser emphasized the 
structure of the output fields. To gain a full understanding of the 
FSF laser it is necessary to evaluate various spectral properties. 
For the field emerging from the FSF cavity the significant observa- 
ble is the optical spectrum. When the FSF laser is used for 
metrology, the significant observable is the RF spectrum from a 
Michelson interferometer. To quantify either of these it is neces- 
sary to have a theoretical description of various correlation 
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functions. Specifically, we require the second-order correlation 
functions for the optical spectrum, considered in this section. We 
require fourth-order functions to describe the RF spectrum, dis- 
cussed in Section 4. 

We begin by considering the FSF laser seeded only by SL This 
provides a simple foundation for more elaborate models of FSF-la- 
ser operation. 

3.7. Passive frequency-independent cavity 

To simplify the problem as much as possible, we first consider a 
passive-cavity model* one in which there is no gain and in which 
losses are independent of frequency. We assume that the cavity 
field originates with the intracavity spontaneous emission. The 
lack of frequency dependence allows usto construct the solution 
using the fact that after one round trip the field E(f) reproduces, 
with no distortion, the field at the earlier time r - x T . The results 
of this simplified approach clarify some of the apparent peculiari- 
ties of the FSF-laser field. 

For our idealized cavity the main losses occur from the small 
nonideality of the AOM, whose efficiency & = exp(-2/ ni ) < 1 is 
taken here to be independent of frequency. We assume that the 
spontaneous emission source in the cavity creates the stochastic 
electric field with the amplitude £ sp (t) and the carrier optical fre- 
quency eo sp . We take £ Jp (f) to have mean value zero and the corre- 
lation behavior 



(22) 



To avoid the unphysical occurrence of an infinite energy and 
infinitely broad spectrum we take the function F(t) to be 



F(t) = !f exp(-r sp |T| + Ha),,, - a>/]T), 



(23) 



rather than the usual delta function F(x) = 6(x). Here r sp and a>sp 
are, respectively, the spectral width and the central frequency of 
the spontaneous-emission line, much larger than all others frequen- 
cies in the system. (For the moment we disregard gain; the assump- 
tion of very large r, p is equivalent to our earlier assumption of 
constant gain within the spectral range of interest.) 

Because we assume here that the cavity losses are independent 
of frequency we can write the amplitude <f (t) of the electric field in 
the cavity as the sum of the immediate spontaneous-emission con- 
tribution t) sp (t) together with emission that was present at the ear- 
lier time r - r r along with similar emissions at a succession of 
earlier times, each diminished by a diffraction & but with un- 
changed time dependance. 

m = 6,(0 + £{ S p(t - r^e-'" + tf^t - 2T r )e- ,A «- T '> M " 

+ &Z sp (t - 3 Tr )e- ,A(f - 2Tr,f -' A(f - Tr)f - ,Ar + . . . (24) 
Thus the cavity field is an infinite sum 



^(f) = s e sp(0 + £^sp(t-nT r )e- 



(25) 



where the phase of component n is <P„ = Jj"o'AT f . 

We infer from Eq. (25) that the spectrum of the field is the sum 
of spontaneous-emission spectra weighted with .ar 2 " and shifted by 
nA. We verify this by examining the second-order correlation func- 
tion for the electric field amplitude, 



C (2) (t, t + t) = <<f (t)<r (r + 1)) exp(i<u/T) 



(26) 



The term exp(io>/T) provides a shift to the optical frequency. 

it is customary to deal with stationary processes for which cor- 
relation functions depend on the time interval r but not on the 
time r. For such processes one can evaluate the power spectrum 



as the Fourier transform of the correlation function; that is the con- 
tent of the Wiener-Khinchine (WK) theorem. For the system con- 
sidered here this theorem does not apply. Instead, we use Eq. 
(25) and the correlation behavior of the spontaneous-emission 
source. Eq. (22), to obtain the result 

C (2, (t, t + t) = i 2 0 £ &~ k F(T -(n - k\r r ) x exp[iS(f . n. k) + ico^t) 

n.k=0 

(27) 

where 



S(f , n ; k) = -(n - k)At + kAx + <f>„ - <P k . 



(28) 



This correlation function C ,2) (t. f + t) depends explicitly on time 
r. This means that the electric field in the FSF laser is not a station- 
ary process. However, this time dependence is periodic, with peri- 
od 2rc/A. Thus we have here an example of a stochastic process 
which is periodically correlated (or cydostationary, periodically sta- 
tionary [30]). For such a process the spectrum can be defined as the 
output of a narrow-bandwidth optical filter. The output of such a 
filter, for such a process, is obtained from the Fourier transform 
of the period-averaged correlation function [30|. When we average 
over time the terms of Eq. (27) with k * n cancel, leading to the 
result 



g (2) (x) = G (2) (t, t + T) = (I £^F(T)e* AT+fa, * r . 



(29) 



k=0 



The laser spectrum is the Fourier transform of this time-aver- 
aged correlation function, 

Ji2){cu) = in 2^ /I dTgi2){T) ex P(- ,C0T )- < 30 > 

This can be considered as an extension of the WK theorem to cydo- 
stationary processes. 

Because the correlation function g {2) (x) is the sum of indepen- 
dent terms, so too is the spectrum the sum of separate spectra. 

#*a«(<o) =ii§if_l dTF(T)e it ° vX ~ ian , (31 ) 

each offset by an increment A and weighted by a power of # 2 , 
/ V) = f> 2 *J^poni(w - a* - JcA). (32) 

k=0 

This result confirms our expectation based upon Eq. (25). 

3.2. Frequency-dependent effective gain 

We next consider a more realistic model of the FSF laser by 
including gain and frequency-dependent losses. The field is given 
by Eq. (19). We assume that the width r sp of the spontaneous- 
emission spectrum is much larger than the width T/ of the intra- 
cavity filter. Thus, the spontaneous-emission process c, p (t) can be 
considered here as a delta-correlated (Wiener-Levy) stochastic 
process: 

{is P (hH; p (t 2 )) = (ls(h-t 2 ). (33) 

When we use this property of c JP (r) to calculate the field, using 
Eq. (19) we obtain the correlation function 



x exp 



x exp 



1=0 

I 



+i(07 + co f )Sx + iMAt + j AM 2 x r - ^ AMi r j . 

(34) 
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We have here defined M and <5t by expressing t as r = Mr r + <5t 
with -T r /2 < <5r < x r /2. 

As in a passive FSF cavity, the correlation function C [2) ((,t + x) 
given by Eq. (34) depends explicitly on time r. As discussed above, 
to obtain the spectrum we need the period-averaged correlation 
function. When we average over time the terms of Eq. (34) with 
M * 0 cancel. The resulting average vanishes unless -r r /2 < r < 
T r /2. when it gives the result 



g (2) (r) = -^r ( X dm £exp(S„ + i(ro + a*)Tl. (35) 
The amplitude of component n here depends on the exponent 

S„ = 2 £>««-/(!&- to)]- (36) 
/-o 

Because the spectral width of J(m) is much larger than 2rc/T, 
this averaged correlation function is just the Fourier transform of 
the spectral density. 



^g (2, (T) = j~ dro/V)exp|+i(ro + a*)t], 



with 



'V)=^£ex P (S n) . 



(37) 



(38) 



This result was obtained earlier (19) using the definition (15) of the 
optical spectral density. 

4. Seeding by spontaneous emission. The radio- frequency 
spectrum 

As we have shown above, the FSF laser output is an example of a 
cyclostationary stochastic process |30], This property has no sub- 
stantial consequences for the second-order correlation functions, 
such as those associated with the optical spectrum, other than 
the need to average over a period. However, the cyclostationarity 
has significant consequences in the fourth-order correlation func- 
tions. These determine, for example, the RF spectrum of the FSF 
output intensity. In this sectipn we analyze such signals. 

4. 1. Passive, frequency-independent cavity 

The output intensity of the idealized FSF cavity, seeded only by 
spontaneous emission, will obviously be noisy. The RF spectrum of 
the fluctuations of the output intensity is of special interest for 
applications to metrology, where it is observed the output of a 
Michelson interferometer. Observations of the RF spectrum were 
important in the use by Nakamura et al. of a FSF laser, not seeded 
from an external source, for ranging measurements [5,6]. RF spec- 
tral measurements continue to be central to the use of externally 
seeded FSF lasers for ranging [9|. 

The RF spectrum derives from the second-order correlation 
function of the intensity, i.e. the fourth-order correlation function 
of the electric field of Eq. (25).The intensity is 

/(0= ^'(0^(0 

= ^ jr # n+k t ip (t - nx r )C 5p (t - kx r ) x eM-war-i*^ (39) 

We define the correlation function for this intensity, based on 
deviation from the mean value 7, as 

G* 4, (r, r) = ([l(t) - 7][/(t + t) - J]) = </(t)/(t + r)> - P. (40) 
We express the averages as the multiple sums 



(l(t)l(t + T)) = (C/8/l) 2 £ 

nJr.mj=0 

* (i ip (t - nx r )(l p (t - kx r )( sp {t + r - mxr)Z' ip 

x (r + T - sx r )} x exp [i(n - k)M + i(m - s)A(f + t) 



(41) 



and 



P = (C/87C) 2 



£ ^ k (is P (t - nx r )(l p (t - to,))*-"*-*-* 



nM=0 



(42) 

To average the product of four random functions we need to in- 
voke additional assumptions about the statistical behavior of the 
function <J jp (t). We will assume that this is a Gaussian random pro- 
cess [31]. This assumption allows the factorization 

(«sp(tl)^(t2)«sp(t3)^(r 4 )) = (fsp(t,){; p (r 2 ))({sp(t3)^(t4)> 

+ (« v (ti)«i(t«))(< v to)«i(ta)). 

(43) 

The width of the spontaneous-emission spectral distribution is 
much larger than all the relevant frequency parameters. Therefore 
we assume that the x dependence of the function F(r) in Eq. (23) is 
much more localized than are the variations with x in Eq. (41 ). We 
carry out the ensemble average of the particular products shown 
above and find that the fourth-order correlation function is inde- 
pendent of time: 



G (4) (t,T) = C (4) (T) = C 0 £ \F(X - MX r )\ 2 & m 

M— no 

with 



(44) 



(45) 



This correlation function does riot depend on the frequency shift 
A. Therefore particular characteristics of the FSF behavior is not 
important; we obtain the same result by considering a cavity that 
has no frequency shifter. The periodic structure of the correlation 
functions occur because the field repeats itself after one round trip. 
The field at any time results builds on superpositions of many con- 
tributions, each shifted in time by x r and weighted by the reflection 
coefficient The interference terms involving different compo- 
nents do not contribute to the mean intensity but they do give 
important contributions to the correlation function. This correla- 
tion function is independent of time and so. by the WK theorem, 
the RF spectrum ] {A) (G) is its Fourier transform. 

f 4) (Q) = ^Jl dxC«\x) exp(-ifli). (46) 

Thus the spectrum of output intensity fluctuations is the infinite 
sum 

f\Q) = g /* dxexp(-iQx) £ \F (T - Mx r )\ 2 ^ M K (47) 

We write this as the product of a scaling factor G 0 and two fre- 
quency-dependent functions, 

f 4) (Q) = C Q L ip (Q)<?(Q). (48) 
Here the smooth envelope 

L ip (Q) = -L /°°dT|F(T)| 2 exp(-ii2T) = 1 fsp/4 2 (49) 
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incorporates the broad but finite width of the spontaneous-emis- 
sion spectrum. For the RF spectrum it suffices to use the approxima- 
tion L ip {Q) = L, P (0) = r sp /Sn. The sum over M in the spectrum (47) 
produces the factor 

*W , t ^expHflMT,) = _ J;f ■ (50) 

M=-3C 

This is a periodic function of Q with a period 2n/x r . For a high 
quality cavity, one with 1 - if 2 <r 1. this function describes a set 
of narrow lines, each centered near a value O q = 2nq/x r for some 
positive integer q. The line profile ¥(Q - Q q ) of each individual 
component is a Lorentzian 

where T c = 2(1 - £)/T r . 

Thus the RF spectrum of the passive cavity output comprises a 
sequence of equidistant spectral features, each with the same nar- 
row profile se(x) t and each centered at an integer multiple of the 
cavity axial-mode frequency-difference 2n/x r . The width of each 
peak is determined by the lifetime of a photon in the cavity. The 
amplitudes of the peaks differ very little. The origin of this spec- 
trum regularity is a periodic repetition of any initial fluctuation 
of the spontaneous emission source. Note, as remarked above, that 
the FSF behavior of the cavity has no effect on these properties of 
the RF spectrum. 

Frames (a) and (d) of Fig. 2 illustrate the properties of the cor- 
relation function Q(x) of Eq. (44) and the spectrum S(Q) of Eq. 
(47). For this figure the various times have been chosen to illus- 
trate the essential physics rather than to describe any actual real- 
ization. The individual peaks of frame (a) have widths set by the 
atomic spontaneous emission lifetime; in practice those widths 
would be much less than the spacing, and the figure would show 
only vertical lines. The dashed-line envelope has a decay time set 
by the lifetime of a photon in the cavity; in practice this decay 
would be much more gradual than is shown. The spectral proper- 
ties of frame (d) are similarly exaggerated. The widths of the indi- 
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vidual peaks are set by the cavity decay time, while the dashed-line 
envelope expresses the spontaneous-emission lifetime. 

4.2. Michelson interferometer with input from a passive cavity 

Here we analyze the RF spectrum of the output of Michelson 
interferometer. For simplicity we consider the interferometer 
formed by an ideal beam splitter and two totally reflecting mirrors. 
We take the difference of interferometer arms to be the distance 
L = cT/2. where T is a total time delay. 

The electric field in the output arm is 

W0«jI'(fl + *<r-T)|. (52) 
The interferometer output intensity, 

l**to-3fi\*(t) + *{t-T)\\ (53) 
can be written as a sum 

fc«(0 = J WO + '(t - T) + Vr(t, t - T)}, (54) 

where l(t) is the output intensity given by Eq. (39). The interference 
term 

/i„ f (t,t-T) = 29?[W("r-r)] (55) 

makes no contribution to the mean intensity but, as we will show, it 
has important effects on the RF spectrum of the Michelson-interfer- 
ometer output. The correlation function for the intensity of Eq. (54), 

T) = (WOWf + T)> - ? Mich (56) 

has two distinct contributions, which we denote as CjJ/ ch (t,T) = 
G^ 4) (t, r) + Cf^t, t). The first of these is the sum of three individual 
correlation functions 

C< 4) (t, T) = G< 4, (t) = ± [2C (4) (T) + C< 4) (T +T) + G (4) (r - T)] , 

(57) 



(5 
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Fig. 2. Left-hand frames (a), (b), (c) show the correlation function, right-hand frames (d). (e), (f) show the resulting spectrum. Top row. frames (a) and (b) show the second- 
order correlation function of Eq.{ 44) and the corresponding optical spectrum of Eq. (47). Middle row, frames (b) and (e) show the correlation function C^ 41 of Eq. (57) and the 
corresponding RF spectrum f? of Eq. (62). Bottom row. frames (c) and (d) show the correlation function O^' of Eq. (58) and the corresponding RF spectrum of Eq. (64). The 
observable spectrum is the sum of contributions (e) and (0- 



306 



LP. Yatsenko et aL/ Optics Communications 282 (2009) 300-309 



and is independent of time. The second part originates with the 
interference intensities, 

CJSce.r) = ~(hm(t,t - T)i ini (t + t,r + 1 - T)). (58) 

This interference term can be averaged in the same way as was 
done when obtaining C i4) (t) in the previous subsection. The result is 

Cfitt. t) = C(T) = Co J F(T - Mx r ) 2 & m cos[MA71 (59) 

The additional phase shift AT between consecutive components in 
the cavity leads to some modulation of the correlation function with 
a period 27rr r /AT. 

The RF spectrum of the output from the interferometer receives 
contributions from two parts, corresponding to the two parts of the 
correlation function. 

/aLw-^w+JJSw- (so) 

The contribution Ji 4) (G) comprises a sequence of equidistant spec- 
tral profiles, each with the same profile <e(Q) but with different 
amplitude: 

x f c/t[2G ,4) (t) + C (4, (t + T) + C* 4) (t - 7)] exp(-iftT) 

x [2F(T - MT r ) 2 + F(T + 7 - Mt r ) 2 + F(T - T - Mr,) 2 ] . 

(61) 

This function is expressible as the product of a scaling factor and 
three frequency-dependent functions {in contrast with the two 
functions of Eq. (48)], 



ft\tot) =^L sp (Q)#(Q)cos 2 (QT/2). 



(62) 



The first factor, L sp (Q), can be considered constant, in keeping 
with our assumption of large r sp . The second factor. ¥(Q), is that 
of Eq. (48): it describes a set of narrow peaks, centered at the dis- 
crete frequencies Q q = 2nq/x r . In contrast to the previous section, 
here there occurs an additional frequency dependence: amplitudes 
of the narrow discrete components are weighted by the factor 
cos 2 (flT/2), evaluated at Q Q . This revision occurs because the cor- 
relation function C 3 originates from an incoherent sum of three 
components delayed by time T. Fig. 2b illustrates the properties 
of the correlation function Gj 4) (t) of Eq. (57) and of the spectrum 
Ji ] {Q) of Eq. (62). The latter exhibits an overall damped cosine 
modulation, as expected. 

The contribution j\*}(Q) to the RF spectrum from Cj^, discussed 
next, is qualitatively different. From the definition we obtain the 
expression 

/i2(0)=*J= £ F(T-MT r ) 2 ^cos(MA7]exp(-ifiT) 

U M=-v 

= ^sp( fl ) £ ^ ,M, cos[MA71exp(rMGT r ). (63) 

M=-oo 

This spectral function therefore comprises a set of doublets. 



Co 



£L sp (Q)<?(Q-2nq/x r -yJ) 



+ §t sp (ft)*(G- 2nq/z r + y c T)]. 



Fig. 2f illustrates this spectral structure. The quantity y c = A/r f 
appearing here is recognizable as the chirp rate used for optical 
ranging. 

To summarize: the total RF spectrum of the Michelson-interfer- 
ometer output consists of narrow discrete components centered at 
the frequencies Q q = 2nq/x ft each weighted by a factor 
cos 2 (fl„T/2). Each of these components is accompanied by a dou- 
blet of components shifted from Q q by ±yj and whose amplitudes 
are almost independent of q. The frequency separation of these 
doublets incorporates information about the difference in length 
cT/2 of the interferometer arms that are used in ranging applica- 
tions of FSF lasers [5-9]. Fig. 2 illustrates, on the third line of 
frames, these features. 

It is worth noting that the existence in the RF spectrum of com- 
ponents shifted by ±y c T is readily understood from a model of a 
FSF laser seeded by a phase fluctuating CW laser [7], Here we see 
that these components appear in a model of the passive FSF cavity 
seeded by uncorrected spontaneous emission as a result of inter- 
ference of cyclostationary stochastic processes. 

4.3. Interferometer input from a frequency-dependent cavity 

Here we consider a more realistic model of the FSF laser by 
including both gain and frequency-dependent losses. The field is 
given by Eq. (19). We will continue with assumptions made in Sec- 
tion 3.2: the width r sp of the spontaneous-emission spectrum is 
much larger than the width r f of the intracavity filter, so we ide- 
alize the spontaneous-emission process £ ip (t) as a delta-correlated 
stochastic process; its correlation function is given by Eq. (33). In 
addition, we assume that the width /> is much larger than both 
the frequency shift A and the cavity axial-mode frequency-differ- 
ence 2lt/X r . 

Using the solution (19) for the FSF laser field and the condition 
r f » A,27r/T r we obtain, after some complicated but straightfor- 
ward algebra, the following expression for the intensity correlation 
function C* 4, (r,T) defined by Eq. (40): 

C' 4) (t,r)=G (4) (r) 

c f x dm 00 2 

= li^ ^ ex ^ iwST) E exp[5(nj,n,M)] . (65) 

Here the argument of the exponent is zero unless 0 < x - Mx r < x r 
when it is 



n-1 



n-1+M 



S(w,n,M) = £>«r -/(t* - /A)] + £ (g Mf -/((Z7 - fA)]. (66) 



1=0 



To simplify this expression, we use the condition T; » A, 27T/r r 
and replace the summation over n in (65) by integration. The expo- 
nent S(nj,n,M) reaches a maximum at 

n = n max - — —21 + M/2, 



with w M = \fml~ M 2 A 2 /4. In the Gaussian approximation justified 

in our earlier paper [19] the parameter m 0 = y/2(g sat -f m )rf is the 
frequency shift between the maximum of the Gaussian fit to the 
optical spectrum J (2) {oj) =J Q exp [-(w -w 0 - o>/) 2 /f 2 ] and the fil- 
ter frequency oj f . As shown in |19], and justified experimentally in 

(20), the shift w 0 and the width f of the FSF-laser spectrum are re- 
lated by the expression 



(64) r 2 ra 0 = y 3 = (Ar 2 ). 



(67) 
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Using the Gaussian approximation we then write 
£ exp|S(ro,n,M)| a Aexp|S(n»,n 1 » 1B ,M)). (68) 



The phase and amplitude occurring on the right hand side are 



3y J 



M 

AMI 



(69) 



A = yfny?. (70) 
Here 

M 0 = f/A » 1 (71) 

is the effective number of round trips of a SE photon in the FSF-laser 
cavity. Integration over m gives finally 

M 2 



C (4, (T) = C, £ F,(T-Mt r ) 2 exp 

M=-oo 



where 
C, 



exp 



8wl 



F ( (T)=^=exp 



r 2 f 2 



(73) 
(74) 



The correlation function (72) is very similar to the correlation 
function of the output intensity of the passive FSF cavity, Eq. 
(44), but there are two differences: 

1. When gains and losses are present the width of the narrow 
function F/(t) is determined by the reciprocal of T , the width 
of the optical spectrum. In the passive cavity this width is fixed 
by the spontaneous-emission correlation time. 

2. The successive maxima of the correlation-function decrease in 
keeping with a Gaussian law exp(-M 2 /Mj) with M 0 T r equal 
to the effective lifetime of a photon in the active FSF cavity. 

It can be shown that the same is true for the correlation func- 
tion of the Michelson-interferometer output excited by the active 
FSF laser: the contribution C^t) is given by Eq. (57) with G* 4) (t) 
defined by Eq. (59) and the interference term 



G!„ 4 f ) (T) = C / £ F|(r - Mx r ) 2 exp 



M=-oo 



M 
*4M 2 0 



cos[MAT]. 



(75) 



Thus we conclude that the RF spectrum of the FSF laser output 
intensity and the Michelson-interferometer output have the same 
behavior as in the passive cavity. The only difference is that the 
amplitudes of narrow discrete components centered at the fre- 
quencies Q q = 2nq/r r and 0^ = 2nq/x T ± y c T decrease according 
to the Gaussian law exp(-(2 2 /r 2 ) and the shape St(Q - Q q ) of each 
narrow component is not Lorentzian but a Gaussian: 



^(x)aexp [-x 2 /(M 0 T r ) 2 ]. 



(76) 



5. FSF laser with seeding by both spontaneous emission and 
continuous-wave radiation 

5.1. The intensity equation 

We consider here an FSF laser seeded simultaneously by spon- 
taneous emission and a monochromatic seed laser. Because there 
is no nonlinearity in our model (only saturation of gain by the total 
intensity) the laser field is the sum of two fields described by Eqs. 



( 1 9) and (20), each created by one of the sources. These sources are 
statistically independent and so Che averaged intensity is also a 
sum of spontaneous and discrete spectrum intensities. 

/ = /sp+/*ed. (77) 

To evaluate the contribution of the external seed, we con- 
sider a monochromatic seed, meaning that its bandwidth is much 
smaller than either the axial mode spacing 2n/x r or the frequency 
shift J. For such a seed the field S it€d (t) is [19] 



<W0 = £, Y^Qn expH<P„ - i(n7 s + nA)t], 



(78) 



n-~0 



where 

* n = -T r n[a> s + (n + l)A/2j. 

We write the constant real-valued amplitudes as 
a„ = exp[S Ifed (g ja|1 n)], where 



(79) 



( 72 ) S^te^.n) 



= n 



f 6w 2 + 6tp,nA + 6n7,A + 2n 2 A 2 + 3A 2 n + A : 



12/7 



(80) 

The contribution l 5Kd of the discrete spectrum to the mean 
intensity / reads 



(81) 



n=0 



where the intensity of the seed laser radiation injected into the FSF 
laser cavity is 



and 



oo 

Fsetdigsat) = exPPWfcanH)]. 



n=0 

The mean intensity of the spontaneous-emission spectrum 
has been evaluated earlier [19]. It can be written as 

hp = Ijpfspfesot)- 

Here 

_cr,j\ 



(82) 
(83) 

I, 

(84) 
(85) 
(86) 



is the intensity of spontaneous emission within the spectral interval 
r f and 



t'/J 



with 

Ssp(gsat> n ) = n 



2{S utt -! m )-~^(n 1 -\) 



- ln(n)/2. 



(87) 



(88) 



Because the stochastic averaging removes any interference con- 
tribution the sum of the two terms (81) and (85) gives the total 
intensity 7, It is this intensity that saturates the gain 

= £o/0 + The determination of / = haiij^; - 1 ) requires 
the solution of the following transcendental equation for g sat : 



'sat 'sor \g,ol J 



(89) 
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Our interest is with FSF lasers in which the output spectrum is 
both smooth and broad. Such conditions are possible if 
'o </«*(*?- 1) and if the spectral filter is sufficiently broad, 
r, » A. For these conditions the general operation of the FSF laser 
is readily understood from Fig. 3. The figure shows the cavity gain 
as a constant, and the loss as varying around the central frequency. 
The frequency shifter displaces any seed frequency toward the 
right. By contrast to the operation of a conventional laser, there ex- 
ists a finite band of frequencies ojf - Wo < (o < (Of + oJo. symmet- 
rically centered about the central frequency of the filter co /( 
where the frequency-dependent effective gain is positive, 
Stax -fm ~ ^ > °- ^ fr e Q uencv component within this band will 
undergo groWth as successive round trips increase the frequency 
by A. However, once the frequency exceeds a>f + ro 0 the losses 
dominate, and the frequency component will diminish as further 
frequency shifts occur. Thus we expect the saturated gain g^ to 
be slightly greater than the minimum Ioss,/ m . 

To proceed we make two approximations. First, because we ex- 
pea the saturated gain g M to be slightly greater than f m we 
approximate it as g^ =/ m on the right hand side of Eq. (89). Sec- 
ond, as in our earlier work (19), we assume that the seed frequency 
is far from the central frequency of the laser spectrum, so that the 
laser output spectrum is very close to a Gaussian form. 

We express S^fg^.n) and S^ig^n) as a Taylor series, cen- 
tered around the value n = n jt where dS,(n)/dr = 0 {i = sp or 
i = seed). and we retain only terms through second order in 
(n - n,). Within this Caussian approximation the number of dis- 
crete components in the interval rxr 0 is large, n7 0 /A» 1, and so 
we replace the summation over n in Eqs. (83) and (87) by an inte- 
gration. We thereby obtain the results 



F^ed as \/nM Q exp 



2w\ + m] - 3n7 5 07§1 



3y3 



and 



exp 



2 w 0 
Eq. (89) now reads 

y 



Ml 

.3)4 



Wq 



exp 



3y 3 J \va 0 



2m\ + ntf - 3m,xol 



3f 



1 

"r 



Here the small parameter p. defined as 



ft-* 



2Ufj-1)y' 



(90) 



(91) 



(92) 



(93) 
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Fig. 3. Qualitative portrait of the magnitude of the FSF laser output-field spectral- 
density as a function of frequency, together with frequency dependence of the gain 
(constant) and loss (quadratic). Two choices of seed frequency are marked. 



describes the spontaneous-emission source. The parameter 

2 r,/ e 



(94) 



expressing the relative importance of the seed laser, is the ratio of 
the intensity of the seed laser to the spontaneous-emission inten- 
sity in the spectral interval y/ny/2. Using experimental data [20] 
we estimate that for an Yb 3+ -fiber laser the parameter e is equal 
to 1 when the power of the seed laser injected into the FSF cavity 
is about 50 uW. 

52. Competition between spontaneous emission and continuous-wave 
seeding 

In this subsection we discuss solutions to Eq. (92). We consider 
two specific choices for the seed frequency. In the first, marked a in 
Fig. 3. the seed frequency is at the minimum of the loss curve. 
m, = 0. In the second, marked b in the figure, the seed is near the 
frequency where growth first becomes possible. 

5.2.1. Resonant seed, a>, = co f 

When the seed frequency coincides with the minimum of 
losses, m % - 0, Eq. (92) takes the form of a quadratic equation 



Y> + eY = - 
involving the variable 

From the solution, 

Y m 

we obtain the intensity of the continuous spectrum 1 $p as 
4//? 



r/»r(>7-l) 



and the intensity of the discrete spectrum as 
2c 



heed : 



'e + y/e 2 + 4/p 



(95) 



(96) 



(97) 



(98) 



(99) 



We see that 7^ > 7 V for e > 1 / >j2p » 1 . This is to be expected: 
when the seed frequency is resonant with the loss minimum, 
o) s = oj ft then the seed laser must compete with spontaneous 
emission that has been amplified throughout the frequency inter- 
val from a)f - Wo to a>/. 

5.2.2. Seed at gain edge, (o s = co f - w 0 

A better choice for the seed laser is o) s = co f - w 0 . Eq. (92) then 
reads 



— exp 
w Q r 



(100) 



The solution to this equation for m 0 with 0 1 gives the following 
expressions for the intensity 1 $p of the continuous spectrum 



7 __L 



l**(ri - 1) 



and the intensity l s „d of the discrete spectrum 



heat = 



1 +ey/d 



W?-i) 



(101) 



(102) 
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where the parameter a is defined as 

This parameter is almost independent of laser parameters, in 
previously reported experimental work [20| the measured value 
of a was a ^ 2. The discrete spectrum dominates in the case of a 
detuned seed laser starting with much smaller intensities: 
>h P fore > 1/v/^~ 1. 

We see that even for optimal detuning of the seed laser there 
exists a large background of continuous radiation originating from 
spontaneous emission. For a strong seed laser, e » 1, the ratio of 
spontaneous to seed intensity is 

j^-^-L (104) 

This means that, roughly speaking, the ratio of continuous to dis- 
crete spectra is equal to ratio of seed laser intensity to the intensity 
of spontaneous emission in the spectral interval w 0 . 

6. Summary and conclusions 

We provide here the first complete theory of the spontaneous- 
emission seeded FSF laser based on correlation functions. We pres- 
ent the second-order correlation function, from which one obtains 
the optical spectrum of the radiation emerging from the FSF cavity. 
Although the correlation function is not stationary, it is periodic, 
and so a period-averaged correlation function leads to a simple 
expression for the spectrum. 

We also present the fourth-order correlation function analysis 
with which one can evaluate the RF spectrum of the laser output 
intensity and of the output of a Michelson interferometer, as is 
used in optical-ranging metrology. We show that even in a FSF la- 
ser seeded solely by spontaneous emission this spectrum com- 
prises a set of doublets, whose spacing gives directly a measure 
of the length difference between the interferometer arms. The exis- 
tence of this structure in the RF spectrum of a the Michelson-inter- 
ferometer output results from correlation of interference terms of 
individual components of a cyclostationary stochastic processes, 
and is not seen in the optical spectrum of the FSF laser. 

We discuss the operation of a FSF laser whose field originates 
from the joint action of spontaneous emission and a continuous- 
wave monochromatic laser seed. We show that in the FSF laser 
the competition between SE and external seed does not act in 
the usual way. of exponential growth of modes with greater gain. 
Instead, there occurs a balance between the two fields that is re- 
lated linearly, rather than exponentially, to the control parameters. 
Although it would be desirable to suppress the spontaneous-emis- 
sion component of the FSF laser output, thereby improving the sig- 
nal-to-noise ratio for optical-ranging technique based on phase 



modulation of the CW seed laser, we show that this is not possible. 
This behavior differs from that of a normal laser. 

We have not here discussed the temporal properties of the laser 
output. The statistical approach we have used in this paper, based 
on an analysis of the correlation functions and the spectral charac- 
teristics, is necessary because we consider a FSF laser seeded by 
spontaneous emission - a sequence of events that is completely 
chaotic. A time series of the laser output, just like the time series 
of the seeding spontaneous emission, is chaotic and conveys no 
information. 
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